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ABSTRACT 

We show that the recently proposed matrix model for M theory obeys 
the cyclic trace assumptions underlying generalized quantum or trace dynam- 
ics. This permits a verification of supersymmetry as an operator calculation, 
and a calculation of the supercharge density algebra by using the general- 
ized Poisson bracket, in a basis- independent manner that makes no reference 
to individual matrix elements. Implications for quantization of the model 
are discussed. Our results are a special case of a general result presented 
elsewhere, that all rigid supersymmetry theories can be extended to give su- 
persymmetric trace dynamics theories, in which the supersymmetry algebra 
is represented by the generalized Poisson bracket of trace supercharges, con- 
structed from fields that form a noncommutative trace class graded operator 
algebra. 
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Recently Banks, Fischler, Shenker, and Susskind [1] have suggested a supersymmetric 
matrix model (based on earlier supersymmetric quantum mechanics models [2]) for use as a 
tool to study uncompactified eleven dimensional M theory, and more recently Banks, Seiberg, 
and Shenker [3] have computed the supersymmetry charge density algebra in this model. 
Although they performed this computation using the Poisson brackets for the individual 
matrix elements, the fact that the final answers are expressible directly in matrix terms 
suggests that it should be possible to perform the calculation using matrix or operator 
methods throughout. We shall show in this paper that the operator dynamics that we have 
proposed [4] and studied with various collaborators [5, 6] is admirably suited for this purpose. 
In the original papers we termed this operator dynamics "generalized quantum dynamics", 
but A. Millard in his thesis [7] uses the briefer and more descriptive name "trace dynamics", 
which we shall use henceforth. What we shall show here is that trace dynamics permits 
an operator treatment of the matrix model for M theory, and generalizes this model to the 
case in which the matrix elements are themselves noncommutative matrices (as would be 
of interest for renormalization group block spin schemes). More generally, we shall show 
elsewhere [8] that all rigid supersymmetry theories can be extended to supersymmetric trace 
dynamics theories, giving supersymmetry representations over trace class noncommutative 
graded operator algebras. Because of the current interest in the model of [1], and the fact 
that this model does not require the gauge fixing technicalities that enter into the case of 
supersymmetric Yang-Mills theory, it is well suited to a brief account illustrative of our 
general approach. 

Let Xi and X2 be two N x N matrices with complex (or more generally, complex 
Grassmann even) matrix elements, and Tr the ordinary matrix trace, which obeys the cyclic 
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property 

TrX 1 X 2 — ^(^l)mn(^2)nm — ^(^2)nm(^l)tim — TrX 2 Xi . (la) 

m,n m,n 

Correspondingly, let 0\ and 6 2 be two N x N matrices with complex Grassmann odd matrix 
elements, which anticommute rather than commute, so that the cyclic property for these 
takes the form 

Tr^x^ 2 = ^(^l)mn(^)nm = ~ ^](^)mn(^l)mri = — Tr# 2 #l • (16) 

m,n m,n 

The cyclic properties of Eqs. (la, lb) are just those assumed for the trace operation TV of 
trace dynamics (although in Refs. [4-6] the fermionic operators are realized as matrices with 
complex matrix elements, all of which anticommute with a grading operator (— 1) F ); we 
shall continue here to use the notation Tr, with the understanding that fermionic matrices 
are Grassmann odd matrices obeying Eq. (lb), while bosonic matrices are Grassmann even 
matrices obeying Eq. (la). From Eqs. (la) and (lb), one immediately derives the trilinear 
cyclic identities 

TtX 1 [X 2 ,X 3 ] =TrX 2 [X 3 ,X 1 ] = TrX 3 [X 1 ,X 2 ] 
TrX ± {X 2 ,X 3 } =TrX 2 {X 3 ,X 1 } = TrX 3 {X u X 2 } 

(lc) 

TrX{9 1 ,e 2 } =Tr6 1 [e 2 ,X] = Tre 2 [9 u X] 
TYM^A} =Tr{e 1 ,X}9 2 = Tr[9 u 9 2 ]X , 
which are repeatedly used below. 

The basic observation of trace dynamics is that given the trace of a polynomial P 
constructed from noncommuting operator variables, one can define a derivative of the number 
TrP with respect to an operator variable O by varying and then cyclically permuting so that 
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in each term the factor SO stands on the right, giving the fundamental definition 



5TiP = Tr 



5TtP 
~50~ 



50 



(2a) 



or in the condensed notation with P = TrP, 




(26) 



In Refs. [4, 5] it is shown that using this definition, one can construct a complete Lagrangian 
and Hamiltonian dynamics for systems with noncommuting graded operator variables and 
a trace Lagrangian, in which a generalized Poisson bracket plays the role played in classical 
mechanics by the classical Poisson bracket, or by the commutator in quantum mechanics. 
In Ref. [6] it is further shown that this dynamical system has the remarkable property that 
its statistical mechanics gives complex quantum field theory, with ensemble averages of the 
operator variables effectively obeying standard canonical commutation relations. 

Let us now turn to the matrix model for M theory. It is described by the trace 
Lagrangian L given by 



with D t O = d t O — i[A ,O]. In Eq. (3a), a summation convention is understood on the 
indices i,j which range from 1 to 9; A and the Xj are self-adjoint N x N bosonic matrices 
with complex number matrix elements, while 6 is a 16-component fermionic spinor each 
element of which is a self- adjoint N x N complex Grassmann matrix, with the transpose T 
acting only on the spinor structure but not on the N x N matrices, so that 9 T is simply the 
16 component row spinor corresponding to the 16 component column spinor 9. The potential 
A Q has no kinetic term and so is a pure gauge degree of freedom. Finally, the ji are a set of 



L = Tr (^DtXiDtX* + i6 T D t 6 + ^[X h X^X*, X>] + Q T li \Q, X 1 ] 



(3) 
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nine 16 x 16 matrices, which are related to the standard 32 x 32 matrices T M as well as to 
the Dirac matrices of spin(8) as described in Danielsson, Ferretti, and Sundborg [9]. This 
finishes the specification of the model. 

Before turning to a study of the model's dynamics and supersymmetries, let us 
summarize some of the properties of the real, symmetric matrices % that are needed. These 
matrices satisfy the anticommutator algebra 

{li,lj} = 2Sij , (4a) 

as well as the cyclic identity (which follows by projection from Eq. (4. A. 6) of Green, Schwarz, 
and Witten [10]), 

£ (S mn 5 pq - -fTlf 9 ) = , (46) 

cycle p-^q^n^p 

with i again summed over and with the indices m, n,p, q spinorial indices ranging from 1 to 
16. Defining 

Hj = \h,lj} , (4c) 

so that 

lilj = Sij + 7^ , (4d) 

one readily derives from Eq. (4b) an identity given in Ref. [3], 

T^Tf + Hh? n + (rn~p) = 2(^S mp - ^8^) . (5a) 

By standard gamma matrix manipulations using Eq. (4a), one also derives the the fact that 
the matrix 

A ijk = liljlk - hjjk + Sik^j - Sjkli (56) 
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is totally antisymmetric in the indices i, j, k (it is just the antisymmetrized product 7[i7j7&] 
with normalization factor |), as well as the identity 

^{lim, lij) = liamlilj] + hj$im - 5 mj S ie , (5c) 

with the first term on the right the antisymmetrized product including normalization 
factor t^j. 

Now let us turn to dynamics. From the trace Lagrangian of Eq. (3), using the 

definition of Eq. (2) we find the operator derivatives 

5L 



SA 



i[X\D t X i }-29 1 9 



5L 

6(dtXi) 



" L - i[D t X\ A ] + [[X^X%Xj] + 29 T f9 



~-D t X l (6a) 



6 -^ = -zD t 9 T + [9 T ,A ]-2[9 T ^X l 
6L -i6 T . 



Substituting these into the operator Euler-Lagrange equations [4] 

d 5L 5L 



dt 8{dtO) SO ' 

and regrouping terms, we get the equations of motion of the matrix model 

D 2 t X l =[[X\ X%Xj] + 29 T f9 

D t e T ^[e^x*] D t 9 = ifr i e,x i ] , 

together with the constraint 



(66) 



(7a) 



C=[X\ D t Xi] - 2i9 T 9 = . (76) 
To transform the dynamics to trace Hamiltonian form, we define the canonical momenta px t 



and pe by 

T' ] T (*») 

so that the trace Hamiltonian is given by 

H = Tr(p Xi d t Xi+p e d t 9) - L = Tr Q OTx * - |[^, X 3 ][X\ X'\ + ip e7i [0, X 4 ] + iA Q C^j . 

(8b) 

Note that because pg = i6 T , it is necessary to write the trace Hamiltonian in a form that is 
manifestly symmetric under the replacements pg — > i9 T , 6 — > — zpj; it is easy to check that 
by virtue of the cyclic identities of Eq.(lc) and the symmetry of 7$ that Eq. (8b) has this 
symmetry. The Hamilton equations following from the trace Hamiltonian of Eq. (8b) are [4] 

dtXi 



dtPx * = ITy 7 

£ 0) 

dm = -le ' 

and it is easy to check that they are the same as the operator Euler-Lagrange equations 
derived above. 

So far we have reproduced standard results of the matrix model, but have broken no 
new ground. Now let us consider the variation of the trace Lagrangian, calculated from 



SL = ^ { ^ SAe + ^ mx , ) + ^ sx , + ^ )Sm + p e] , (10o) 



for the supersymmetry transformation defined by 

5X l = - 2e T 1 l 9 = 29 T fe 

59 = - ^S + ^[^^ J h^ + e' (106) 
5A = - 2e T 9 = 26 T e . 
Here e and e' are 16 component Grassmann c-number spinors, that is, they are column 
vectors each of whose 16 components is an independent lxl Grassmann matrix. Using 
only the cyclic trace identities and gamma matrix properties given above, it is a matter 
of straightforward but lengthy calculation to verify that the trace Lagrangian is invariant 
under the transformation of Eq. (10b) when e and e' are time independent. Note that in 
this calculation the variables X iy 9, and A are treated simply as noncommuting operators 
that are unspecified apart from their bosonic or fermionic character; in particular, we do not 
have to assume that their matrix elements when they are written as N x N matrices are 
c-numbers or Grassmann c-numbers. In other words, we have shown that the matrix model 
of Eq. (3) is still supersymmetric when viewed as a model over general noncommutative trace 
class graded operator variables. 

When e has a time dependence 5L is no longer zero, but instead is given by 

5L =d t Tr[- i9 T e' + {9 T li D t X i + \i9 T ll0 [X\ X j ])e] 

2 (11a) 

+Ti[2i9 T d t e> + {29 T li D t X i - ie T lij [X\X^])d t e] . 
This identifies the trace supercharges Q^, and Q a as 

=Tr2i9 T a 

(116) 

Q a =Ti(29 T ll D t X t - %9 T ll3 \X\ X'])a , 
and their conservation is easily checked using the equations of motion and constraint of 
Eqs. (8a, b) and the identities of Eqs. (4, 5). To check the supersymmetry algebra, we must 

9 



first write the supercharges of Eq. (lib) in Hamiltonian form, symmetrized with respect to 
po and iO T , giving 

=Tr(p + i6 T )a 

1 (12«) 



Q a = - Tr(p e + i6 T ) (i liPXi + - ltJ [X\ Xi}) 



a 



In addition to the supercharge Q Q , Ref. [3] also introduces a supercharge density; in our 
language this can be written in terms of a Grassmann spinor operator (3 as 



Q^ = -^r{(p0 + z9 T ),t lt p Xi + ^ lt3 [X\X^}(3 . (126) 



Let us now introduce the generalized Poisson bracket [4, 5] corresponding to the 
Hamiltonian structure of our model, defined for any A = TrA and B = TiB by 



{A,B} = Tr( 



SA 5B 5B SA 5A5B SB5A\ 

W i S P ^ i ~dX l 6p^~l9 6p~e + l9Sp~e) ' (13) 



We can now use the generalized Poisson bracket to give a basis-independent evaluation of 
the supersymmetry charge algebra. For the brackets involving Q', we easily find 

{Q^,Q^} = - 2tTra(3 

1 (14) 
{Q' a , CM =2?Tr(iT?>p Xi + -1*[X\ Xi])a a p b . 

The only case involving significant work is the bracket of Q a with Q^, which on 
substituting the operator derivatives of the Q's into Eq. (13), but before further algebraic 
rearrangement, takes the form 

{Qa,Q/3} = Trf/ , (15a) 
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with U given by 

U=2[Xt,a b 2W a ±1$][p d ,i9 c ]vy? d 

+[X m ,[2i9 c ,p d ]]^a b 2ie a iryf 

1 1 1 (156) 

+ \m,ilfpx e + l<[X e ,X m ]}( i7 fp Xi + \^[X\Xi])a b . 
Algebraic rearrangement of this using the cyclic and gamma matrix identities gives TrU = 

TrT, 

T = T x ±ib M a h {3 d + T|2 7 f a 6 ft + X^X^^^^rc^ , (16a) 

with 

2 4 2 (166) 

TrTi =H 

and 

Ti={ Px ,,[x\x^}-z[e a ,[x^,e a ]] 

(16c) 

= - {r,c} + [x\ { Pxi ,xi}\ - i{e a , {x*, e a }} , 

where in getting the final line of Eq. (16c) we have used the Jacobi identities for mixed 
commutators and anticommutators. Equations (14) and (16a-c) agree with the charge density 
algebra computed in Ref. [3]. In the case when (3 is a c-number, we see that the generalized 
Poisson bracket of the trace supercharges gives the trace Hamiltonian, showing that we 
have constructed a supersymmetry representation over noncommutative trace class operator 
dynamical variables. As noted earlier, we will show elsewhere [8] that this construction 
works for general rigid supersymmetric theories, and in particular for the Wess-Zumino and 
supersymmetric Yang-Mills theories in four dimensions. 

We close with a remark on the quantization of the matrix model just described. 
There are two possible points of view one could take. The first possibility would be to 
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quantize by replacing the Poisson brackets for matrix elements used in the computations of 
Ref. [3] by a commutator/anticommutator algebra; in this case one would be treating the 
matrix elements themselves as quantum operators. An alternative possibility is discussed in 
Ref. [6], which is to regard the trace dynamics as fundamental without quantization, and to 
consider its statistical mechanics assuming ergodicity (an assumption that may presuppose 
taking the large N limit; we also remark that the proof of the generalized Liouville theorem 
in Ref. [6] extends, with minor modifications, to the case 

employed here in which the fermions are realized with Grassmann matrices.) It 
is then shown that the statistical averages of the dynamical variables obey the rules of 
complex quantum field theory, with the effective Planck constant given by the expectation 
of the conserved operator C discovered by Millard [7], which for a generic trace dynamics 
model has the form 

C= ^2h,Pi\- iVi'Pi} ■ ( 17a ) 

bosons fermions 

Since for the trace Lagrangian of Eq. (3) C is just the expression of Eq. (7b), which vanishes 
as a constraint, the analysis of Ref. [6] implies that the statistical averages will obey classical 
mechanics, with the X { all commuting with one another and with the corresponding p Xj - As 
discussed in Ref. [4], it is easy to modify the model so that the constraint of Eq. (7b) reads 
instead C = ih; this is done by adding to L the term 

AL = -TrA , (17b) 

which gives a trace action AS = f dtL that is invariant under gauge transformations that 
vanish (or are periodic) at t — ±oo. As we shall show in Ref. [8], it is easy to construct 
supersymmetric trace dynamics theories in which the operator C, although conserved, is 
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not constrained to vanish; the trace dynamics extension of the Wess-Zumino model is an 
example. 
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